Diffusion limited reactions in confined environments 
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We study the effect of confinement on diffusion Hmited bi-molecular reactions within a lattice 
model where a small number of reactants diffuse amongst a much larger number of inert particles. 
When the number of inert particles is held constant the rate of the reaction is slow for small reaction 
volumes due to limited mobility from crowding, and for large reaction volumes due to the reduced 
concentration of the reactants. The reaction rate proceeds fastest at an intermediate confinement 
corresponding to volume fraction near 1/2 and 1/3 in two and three dimensions, respectively. We 
generalize the model to off-lattice systems with hydrodynamic coupling and predict that the optimal 
reaction rate for mono-disperse colloidal systems occurs when the volume fraction is ~ 0.18. Finally, 
we discuss the application of our model to bi-molecular reactions inside cells as well as the dynamics 
of confined polymers. 

PACS numbers: 



It is a somewhat surprising fact that the total concen- 
tration of protein within a cell rivals that within a pro- 
tein crystal[l|. This highly crowded environment plays an 
important role in dynamical processes, such as rates of 
chemical reactions, and thermodynamic properties, such 
as chemical equilibria, observed in vivo]^. While the 
volume fraction of macromolecules within the cell may 
exceed 30%, there may be just a few copies of a given pro- 
tein corresponding to a concentration of a few nanomo- 
lars for a cell of volume l/im"^[J, [^]. This is in stark 
contrast to in vitro experiments where the reactants are 
present at a relatively high concentration with a neg- 
ligible level of crowding molecules. Therefore, caution 
is required when interpreting in vivo biochemical experi- 
ments as it is not always obvious what effect the presence 
of crowding molecules will have. For example, the rate 
of a reaction may be increased if the crowding favors a 
compact transition state, decreased if the reaction is dif- 
fusion limited, or unaffected if the reactants are small 
compared to the crowding species [g. 

In this letter we study the effect of crowding on re- 
action rates in a finite system with a fixed number of 
particles. A biological cell represents such a system. We 
find that the presence of non-reactive particles leads to a 
non-nionotonic reaction rate as the volume of the system 
is changed. This non-monotonicity is the result of two 
distinct dynamical regimes. When the system volume is 
very large the effect of the crowding particles is negli- 
gible. Therefore, the rate of reaction is dictated by the 
time required for the reactants to diffuse the mean sep- 
aration between reactive particles which increases with 
the system size. We will call such systems "concentra- 
tion limited" . At the opposite extreme is the case where 
the system volume is very close to the sum of the total 
volume of the reactant and crowding particles. In this 
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case the separation between reactants may be quite small, 
but the reaction proceeds slowly because the high den- 
sity of crowding particles severely impedes the diffusion 
of the reactants. Such systems are "crowding limited" . 
Between these two limits there is an optimal volume at 
which the reaction proceeds the fastest. 

The competition between the concentration lim- 
ited and crowding limited regimes may be understood 
through the following simple argument. The flux of par- 
ticles at the surface of an absorbing sphere is AiraCrD, 
where a is the radius of the sphere, and the concen- 
tration of particles far from the sphere, c^, is inversely 
proportional to the volume of the system Cr oc i?~^[7|. 
The effect of crowding particles that do not react with 
the absorbing sphere may be included through a rescal- 
ing of the diffusion constant D provided the distance the 
reactant particles must travel is large compared to the 
mean spacing between crowding particles. For the case 
of a lattice-based system, a mean-field diffusion constant 
(which is exact when the crowding particles move much 
faster than the reactants) is achieved by multiplying the 
"bare" diffusion constant Dq by the success rate for a 
particle move. This results in a diffusion constant of the 
form D ^ Dqp where p = (1 — c) is the probability that 
the neighboring site for an attempted move is unoccu- 
pied, and c is the number density of particles on the lat- 
tice. Therefore, the flux at the absorbing surface scales 
as a(l — Nb^R~^)/R^, where N is the total number of 
particles in the system and b is the lattice constant. This 
expression has a maximum when the system size is such 
that the particle density is 1/3. 

In order to go beyond this simple argument, we have 
explored the transition from concentration limited reac- 
tions to crowding limited reactions using an analytically 
tractable model with two reactant particles that react 
instantaneously upon contact. With this simpliflcation, 
the reaction "rate" is just the inverse mean first passage 
time for the particles to find each other. If we make 
the further simplification of holding one of the reactants 



fixed, the mean first passage time is given by 
1 
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where the integral is over the d-dimcnsional volume of the 
system. t{x), the average time for the mobile reactant 
at X to reach the stationary target, satisfies the equation 



Dy^T{x) = -1. 



(2) 



Eq. [5] is subject to a reflecting boundary condition at the 
system periphery and an absorbing boundary condition 
at the surface of the stationary reactant. 

If we place the fixed reactant at the center of a spherical 
box, as shown in Fig. [U Eq. [T] is exactly solvable with 
the result (in two and three dimensions) 
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Here a is the radius of the stationary target, and D{c) is 
the concentration dependent diffusion constant. In the 
spirit of the mean-field argument introduced earlier, we 
employ a lattice model and utilize an excellent approx- 
imation for the self-diffusion of a lattice gas derived by 
van Beijeren and Kutner [lOJ ]. For reactant and crowd- 
ing particles with equal mobilities on a square lattice, the 
self-diffusion constant takes the form 

Die) = ^ (v/4(l - c) + c2(^ - 1)2 - dn 1)) , (5) 

where b is the lattice spacing, F is the attempt rate for 
particle moves, and c is the ratio between the total num- 
ber of mobile particles to the number of accessible sites 
c ~ b^{N + l)/(7r(i?2 - a^)). iFmm Eqs. [3] and E we 
find an expression for the reaction time which is plotted 
in Fig. H 

We have preformed Monte Carlo simulations of this 
model consisting of one reactant and N inert particles 
confined to a two dimensional circular box of radius R 
with a circular target of radius a placed at the center, 
as shown in Fig. [TJ Each run of the simulation begins 
with a random configuration of the reactant and crowding 
particles and ends when the reactant reaches the center 
target at which point the first passage time is recorded. 
The average reaction time is determined from 10,000 such 
runs. 

For the purposes of our simulation we define the box 
and the target to include all sites on the lattice with a 
distance to the origin less than or equal to R and a, re- 
spectively. Although the definition of a circular box on a 
square lattice is somewhat cumbersome, this geometry is 
advantageous relative to square or rectangular boxes be- 
cause the circular box allows for non-integer adjustments 
of the linear box dimension with minimal change to the 



FIG. 1: Snapshot from the simulation described in the text. 
The stationary target is indicated by the black circle in the 
center, and the mobile reactant is the black dot halfway be- 
tween the target and the top of the outer circle. Also shown 
are the 2000 inert crowding particles (grey). 



overall geometry. This fine-tuning of the box size greatly 
increases the number of data points that can be collected 
in the crowding limited regime. 

The simulation used N = 2000 crowding particles 
which, like the single mobile reactant, each occupied a 
single lattice site. The central target had a larger size 
a = 3 to minimize the discreetness effect of the lattice 
when comparing to the continuum theory given by Eq. [T] 
The radius of the confining box was varied from i? = 150 
to R = 26.7 lattice spacings, corresponding to densities 
from c = 0.03 to c = 0.90. The results are shown in Fig. 

m 

The simulation results show a minimum in the reac- 
tion time very near the minimum of 36.9 predicted by 
Eq. [3] (corresponding to 47% of the sites being occupied). 
The reaction times diverge sharply when the box radius 
becomes less than 29 lattice sites, or the concentration 
exceeds ~ 75%. At large box radii the reaction times 
increase as R^ and approach the reaction rate in the ab- 
sence of crowding particles (green line) . As shown in Fig. 
[2] the difference between the crowded and uncrowded re- 
action times is nonzero at all system sizes. This can 
be explained by noting that while the probability that a 
given site is occupied scales as R~'^, the required number 
of steps scales like R"^. The result is that the number of 
time steps where the particle is unable to move due to 
crowding remains nearly constant as the system size is 
increased. 
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FIG. 2: (color online) Comparison of the reaction time pre- 
dicted by Eq. [l] (black dashes) to the simulation data (red 
dots). The green dashes indicate the reaction time predicted 
in the absence of crowding particles. Error bars (not shown) 
are smaller than the size of the dots. 



The reaction times in Fig. [2] deviate from the theoreti- 
cal curve at system sizes below 506 with a maximum error 
of 39% at the highest concentration simulated. This dis- 
crepancy is due to a non-diffusive, logarithmic correction 
to the mean squared displacement (MSD) in the lattice 
gas system [ll|. This correction allows the reactant to 
sample its immediate surroundings more efficiently than 
a purely diffusive particle. The effect of the correction 
term can be seen in Fig. [3] where we compare the MSD 
of the reactant in our simulations to the expected MSD 
(as calculated from the Green's function) for a random 
walker in a circle with diffusion constant given by Eq. [S] 
At the largest system size, when the particle density is 
0.1, the agreement between the simulation and calculated 
MSD is better than 3% at aU times, kl R = 35.8 and 
R = 26.7, where the particle density is 0.5 and 0.1 respec- 
tively, the discrepancy between the observed and calcu- 
lated MSD is 8% and 18% respectively at short times and 
shrinks at later times as the MSD saturates. 

Our results can be extended to the more physically 
relevant three dimensional, off-lattice system with a few 
modifications. We recall that hard sphere systems un- 
dergo a glass transition at a volume fraction (\)c ~ 0.58, 
which is less than the close packing density [9]. Al- 
though the system is not completely frozen, the rele- 
vant time scales diverge sharply. Therefore, we restrict 
our analysis to densities below 0c- In the fluid regime, 
for < 0c; the self diffusion constant is modified not 
only by the short range excluded volume interactions, 
but also by long range hydrodynamic coupling mediated 
by the solvent |l3|. This latter effect is not included in 
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FIG. 3: (color online) The mean mean squared displacement 
of the mobile reactant measured in our simulation (solid lines) 
compared to the calculated mean squared displacement for a 
random walker in a circle(dashed lines). 



our lattice model due to the lack of solvent. Using the 
self-diffusion constant derived in reference [191 (equation 
4.34), together with Eq. |4l we predict that the opti- 
mal reaction time for mono-disperse, hard-core, Brown- 
ian spheres occurs when the system size is such that the 
volume fraction is ^ 0.18. 

Experimentally, a crowded reaction with adjustable 
volume could be realized using microfluidic techniques. 
In this case the tracer and crowding particles would be 
confined to a micro droplet reaction vessel whose volume 
could be controlled through osmotic gradients across a 
semi-permeable barrier [12|. This system mimics in vivo 
experiments that show increased tracer particle diffusion 
in osmotically swollen cells |13l[l4|. 

In most systems of interest the target reactant would 
be mobile rather than stationary. One expects that the 
ratio of the reaction time for the stationary target to the 
reaction time for two mobile reactants in a fixed box size 
should be independent of the number of crowding par- 
ticles present, as the crowding particles merely rescale 
the diffusion time. Our simulations support this intu- 
ition, and we find that the ratio is ^ 1.6 and ^ 2.0 in 
two and three dimensions respectively, with corrections 
on the order of a/ R. 

In the context of these theoretical results it is natu- 
ral to ask to what extent do cells optimize biochemical 
reaction rates by adjusting their size? At first glance it 
would seem unlikely that cell volume would be a useful 
parameter for the cell to use to regulate reaction rates 
due to the large number of reactions that occur simul- 
taneously. Furthermore, differential protein expression 
rates during the cell cycle could modify the overall pro- 
tein concentration and thereby alter the cytoplasmic dif- 



fusion rate[15[. However, it is well known that the cell is 
more complicated than a "bag of enzymes," and there- 
fore, the potential exists for the cell to compartmentalize 
reactions in such a way that the volumes of the compart- 
ments are individually tunable. For example, the diges- 
tion of a pathogen within a macrophage occurs within a 
vesicle created by endocytosis. The volume of this vesicle 
is, in principal, adjustable by the amount of membrane 
used during vesicle creation. Similarly, other membrane- 
bound organelles such as golgi, endoplasmic reticulum, 
mitochondria, and the nucleus could be individually ad- 
justed to optimize reactions occurring within. This is 
consistent with the finding that diffusion rates in the mi- 
tochondria and endoplasmic reticulum can differ substan- 
tially from the cvtoplasmfla [l7| . 

The case of the cell nucleus deserves special consid- 
eration. Here the primary reactions of importance in- 
volve the manipulation of the genome which is encoded 
by DNA. In analogy to our two reactant model, many of 
the genome management functions the cell performs re- 
quire two specific portions of the DNA to find each other. 
These internal cyclization reactions may occur between 
monomers separated by polymer spacers ranging from 
less than a persistence len gth , up to lengths on the order 
of the chromosome size [20, l2l| . If we identify the polymer 
segments that flank and bridge the reacting monomers as 
crowding particles, then we can immediately generalize 
the preceding argument for the rate of crowded reactions 
to the internal cyclization of a confined polymer. Specif- 
ically, the rate of internal cyclization will have a non- 
monotonic dependence on the size of the box containing 
the polymer. This non-monotonicity has been previously 



observed in computer simulations [22l|. 

The onset of crowding limited dynamics at small sys- 
tem sizes leads to non-monotonic behavior in the relax- 
ation of other structural properties of the polymer such 
as the end-to-end vector |23'|. This non-monotonicity is 
quantitatively different from that of crowded reactions 
due to the sub-diffusive behavior of the monomers im- 
posed by the connectivity constraint. However, it can be 
easily explained by noting that the initial effect of reduc- 
ing the volume accessible to the polymer to smaller than 
its unconfined size is to reduce the conformational phase 
space that the polymer can sample. This allows the poly- 
mer to sample the phase space faster resulting in shorter 
relaxation times. However, like the free particle case, as 
the monomer density approaches the close-packing limit 
the polymer becomes jammed and the relaxation times 
lengthen. 

We have shown that a simple, analytically tractable 
model is able to quantitatively predict reaction times 
occurring within crowded environments. This model 
can be generalized to off-lattice systems and systems 
with explicit solvent with slight modification. The non- 
monotonicity in the reaction times shown here has broad 
implications for reactions within cells as well as the dy- 
namics of confined polymers. 

Acknowledgments 



We would like to thank B. Chakraborty, S. Redner, 
and A. Grosberg for useful discussions. This work is sup- 
ported by the NSF Grant DMR-0403997. 



[1] A.B. Fulton, Cell 30, 345-347 (1982). 

[2] S.B. Zimmemian and A. P. Minton, Annu. Rev. Biophys. [15] 

Biomol. Struct. 22, 27 (1993). 
[3] A.R Minton, J. Biol. Chem.276, 10577 (2001). [16] 

[4] S.B. Zimmerman and S.O. Trach, J. Mol. Biol. 222, 599- 

620 (1991). [17] 

[5] D. Hall and A. P. Minton, Biochem. Biophys Acta 1649, 

127-139 (2003). [18] 

[6] A.P. Minton, J. Cell Sci. 119, 2863-2869 (2006). 
[7] H.C. Berg Random Walks in Biology, Princeton Univer- [19] 

sity Press (1983). 
[8] S. Redner A Guide to First-Passage Processes, Cam- [20] 

bridge University Press (2001). 
[9] E. R. Weeks, J. C. Crocker, A. C. Levitt, Andrew [21] 
Schofeld, and D. A. Weitz, Science 287, 287-631 (2000). [22] 

[10] H. Van Beijeren and R. Kutner, Phys. Rev. Lett. 55, 

238-241 (1985). 
[11] M.J. Saxton, Biophys. J. 52, 989-997 (1987). 
[12] J.-u. Shim, G. Cristobal, D.R. Link, T. Thorsen, Y. Jia, 
K. Piattelli, S. Fraden, J. Am. Chem. Soc. 129,8825-8835 
(2007). 
[13] H.P. Kao, J.R. Abney and A.S. Verkman J. Cell Biol. [23] 

120, 175-184 (1993). 
[14] R.Swaminathan, C.P. Hoang, and A.S. Verkman, Bio- 



phys. J. 72, 1900-1907 (1997). 

M.B. Elowitz, M.G. Surette, P. Wolf, J.B. Stock and S. 

Leibler, J. Bacter. 181, 197-203 (1999). 

A. Partikian, B. Olveczky, R. Swaminathan, Y. Li, and 

A.S. Verkman, J. CeU Biol. 140, 821-829 (1998). 

M.J. Dayel, E.F.Y. Hom, and A.S. Verkman, Biophys. J. 

76, 2843-2851 (1999). 

A. Kasper, E. Bartsch, and H. Sillescu, Langmuir 14, 
5004-5010 (1998). 

M. Tokuyama and I. Oppenheim, Physica A 216, 85-119 
(1995). 

B. Miiller-Hill, The Lac Operon: A Short History of a 
Genetic Paradigm, Walter de Gruyter: Berlin, 1996. 
J.E. Haber, Annu. Rev. Genet. 32, 561-599 (1998). 

C. F. Abrams, N.-K. Lee, and A. Johner, Macro- 
molecules, 39, 3655 -3663 (2006). The non-monotonic 
behaviour discussed in this reference is not related to 
that presented in the current work. However, Figure 2 
shows an increase in the internal cyclization rate at small 
system sizes that the authors attribute to "increased re- 
striction of monomer mobility" . 

J. Kalb, J.D. Schmit, and B. Chakraborty in preparation 



